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We analyze the formation and dynamics of bright unstaggered solitons in the disk-shaped dipolar 
Bose-Einstein condensate, which features the interplay of contact (collisional) and long-range dipole- 
dipole (DD) interactions between atoms. The condensate is assumed to be trapped in a strong 
optical-lattice potential in the disk's plane, hence it may be approximated by a two-dimensional 
(2D) discrete model, which includes the on-site nonlinearity and cubic long-range (DD) interactions 
between sites of the lattice. We consider two such models, that differ by the form of the on-site 
nonlinearity, represented by the usual cubic term, or more accurate nonpolynomial one, derived from 
the underlying 3D Gross-Pitaevskii equation. Similar results are obtained for both models. The 
analysis is focused on effects of the DD interaction on fundamental localized modes in the lattice (2D 
discrete solitons) . The repulsive isotropic DD nonlinearity extends the existence and stability regions 
of the fundamental solitons. New families of on-site, inter-site and hybrid solitons, built on top of 
a finite background, are found as a result of the interplay of the isotropic repulsive DD interaction 
and attractive contact nonlinearity. By themselves, these solutions are unstable, but they evolve 
into robust breathers which exist on an oscillating background. In the presence of the repulsive 
contact interactions, fundamental localized modes exist if the DD interaction (attractive isotropic 
or anisotropic) is strong enough. They are stable in narrow regions close to the anticontinuum limit, 
while unstable solitons evolve into breathers. In the latter case, the presence of the background is 
immaterial. 

PACS numbers: 03.75.Lm; 05.45.Yv 



Disk-shaped two-dimensional (2D) Bose-Einstein condensates (BECs) have been created using a superposition of a 
tight optical trap, formed by a pair of parallel strongly repelling (blue-detuned) light sheets, and a loose in-plane radial 
magnetic trap Such configurations may exist with both self-repulsive and self-attractive collisional nonlinearity, 
which can be controlled by means of the Feshbach-resonance technique. In the low-density limit, the dynamics of the 
2D condensate is modeled by the 2D Gross-Pitaevskii equation (GPE) with the cubic nonlinearity. If the density is 
not very small, a consistent reduction of the underlying three-dimensional cubic GPE to two dimensions leads to the 
2D nonpolynomial nonlinear Schrodinger equation (NPSE) 0, Q (a related but different approach leads to another 
form of the NPSEs for ID and 2D configurations, in the case of the self-repulsive nonlinearity [4(). 

A strong 2D periodic optical-lattice (OL) potential, applied to the disk-shaped BEC, splits it into a planar array of 
droplets coupled by weak nearest-neighbor linear interactions, due to the tunneling of atoms across potential barriers 
which separate the droplets [5|. In this situation, the BEC can be modeled by the 2D discrete GPE, derived from its 
continuous counterpart by means of an expansion based on localized Wannier functions @, 0] ■ The same discretization 
procedure, if applied to the continuous equation with the nonpolynomial nonlinearity, leads to the discrete variant 
of the NPSE, as recently demonstrated in detail in the ID geometry A noteworthy feature of the NPSE models, 
both continuous and discrete ones, is that, in the case of the attractive local nonlinearity, they account for the onset 
of the collapse even in the framework of the ID description (recall that the GPE with the self-focusing cubic term 
does not give rise to the collapse in ID). 

The BEC dynamics is strongly affected by long-range interactions in dipolar condensates, which can be formed by 
chromium atoms in an external magnetic field, as shown in experiments @. In particular, the dipole-dipole (DD) 
attraction in the condensate gives rise to a specific (<i-wave) mode of the collapse in the condensate Neverthe- 
less, the dipolar condensate can be stabilized against the collapse, adjusting the scattering length of the collisional 
interactions via the Feshbach-resonance technique [III- It was also proposed to create a condensate dominate d by 
DD interactions between electric dipole moments, that may be induced in atoms by an external electric field [12|. 




A similar situation may be expected in BEC formed by dipolar molecules jl 31 ] - In particular, the creation of LiCs 
dipolar molecules in an ultracold gas mixture has been reported (l4| . 

Once dipolar condensates are available to the experiment, a natural issue is how the nonlocal DD interactions 
affects localized modes (solitons). Various possibilities to create 2D solitons in dipolar condensates were explored 
theoretically. It was demonstrated that isotropic solitons [l5[ and vortices [l6| may exist in the 2D disk-shaped 
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configuration with the moments polarized perpendicular to the disk's plane, if the sign of the DD interaction is 
inverted by means of the rapid rotation of dipoles [ItJ ■ The ordinary DD interaction (with the uninverted sign) can 
give rise to stable anisotropic 2D solitons, if the dipoles are polarized in the plane of the confining disk [l8j]. It may 
be relevant to mention that stable 2D solitons supported by nonlocal interactions (different form those corresponding 
to the DD forces) are known in optical media with the thermal nonlinearity (l9| . 

Solitons in ID models of the dipolar BEC have been recently predicted too, using the respective GPE [20] or NPSE 
[2l| . In the latter case, the model enables to study the influence of the long-range DD nonlinearity on the onset of the 
collapse in the condensate with the self-attractive nonlinearity. Stable ID solitons were also predicted in the model of 
the Tonks-Girardeau (TG) gas of bosons carrying permanent moments, using the known model of the TG gas with 
the self-repulsive quintic local term, to which the long-range DD attraction was added [22| . 

A natural extension of the theoretical studies of localized modes in the dipolar BEC is to include a strong OL 
potential, which leads to the discrete model with the DD interactions between remote lattice sites. Recently, this 
analysis was performed in the ID setting [H, [U, [25[ . In particular, unstaggered solitons in such models, with both 
the cubic and nonpolynomial on-site nonlinearities, were studied, respectively, in Refs. [23| and with a conclusion 
that the DD interactions might enhance the stability of the discrete solitons. 

The objective of the present work is to introduce discrete models for the dipolar condensate trapped in a deep 2D 
lattice, and analyze the existence and stability of localized modes in these settings (discrete bright solitons) . We will 
focus on the most fundamental case of unstaggered solitons, performing the analysis in terms of two different lattice 
models. The first model includes the ordinary cubic on-site nonlinearity, i.e., it corresponds to the 2D version of 
the cubic discrete GPE (to be abbreviated as DGPE), while the second model deals with the nonpolynomial local 
nonlinearity, which corresponds to the 2D discrete NPSE (alias DNPSE). Both models incorporate the cubic nonlocal 
term accounting for the DD interactions. The results obtained for discrete solitons in the frameworks of both models 
will be compared. 

The paper is structured as follows. The two discrete models are formulated in Section II. Families of on-site, inter- 
site and hybrid 2D solitons in the case of the attractive contact interaction, combined with isotropic and anisotropic 
DD interactions, are presented in Section III. Their stability is analyzed in the same section. Properties of fundamental 
localized structures that may be formed as a result of the competition between the repulsive local nonlinearity and 
attractive DD isotropic interactions (i.e., with the inverted sign of the DD force, as mentioned above), as well as 
between the repulsive local and anisotropic DD interactions with the natural sign, are considered in Section IV. The 
paper is concluded by Section V. 



II. THE MODELS 



The scaled form of the continuous GPE, which includes both contact and DD interactions in the 3D geometry, is 
well known 0, : 
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where U(r) is the external potential, 7 = 2Na s /a z is the normalized strength of the local interactions, N the number 
of bosonic atoms, a s the scattering length of inter-atomic collisions, and a z = \Jhj (mw z ) is the scale of the trapping 
provided by potential (1/2) mu> 2 z 2 , which acts perpendicular to plane {x, y) of the disk-shaped condensate. The total 
potential, U(r), includes both the latter term and the in-plane 2D OL potential, W(x,y). The notation used in Eq. 
([T]) assumes that wave function ip (r, t) is subject to the usual normalization, 



\ip{r,t)\ 2 dr = 1. 



(2) 



Coefficient g in Eq. ([T]) defines the strength and sign of the DD interactions between atomic dipoles, and 9 is the 
angle between vector (r — r') and the orientation of dipoles, which is fixed by a strong external magnetic field. 
The reduction of Eq. ^ to the 2D form is performed by assuming the factorization of the 3D wave function 0, S] : 

V'M) = —^=<l>(x,y,t)f{z,r)), (3) 
V 27r M 

where rj is the width of the distribution along axis z (77 may depend on coordinates x, y and time t), <f>(x,y,t) is an 
arbitrary in-plane (2D) wave function, and / (z; rf) is its normalized axial counterpart, 
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With regard to the underlying 3D normalization condition ([2]), the factorized ansatz based on Eqs. Q and Q yields 
the following normalization condition for the 2D wave function, 



\<j>(x,y)\ 2 dxdy = v 7271 " ItI 



(5) 



Substituting the factorized expression (|3|) in Eq. ([I]) , assuming the presence of the strong 2D optical lattice potential 
in the plane of (x,y), and averaging the result along coordinate z, the following 2D NPSE with the nonlocal term 
accounting for the DD interactions is derived: 
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where H is the sign of the contact interaction (H = —1 and +1 correspond for the self-attraction and repulsion, 
respectively), and T = g/ {ypZr; I7Q is the relative strength of the nonlocal (DD) interaction versus its local counterpart 
which can take positive and negative values. In particular, the nonlocal term in Eq. ([6]) is obtained by the integration 
in the direction of z, assuming the long-distance limit (26l.[27j. i.e., that characteristic length scales in the (x, y) plane 
are essentially larger than the above-mentioned transverse-confinement size, a z . 

In the presence of the strong in-plane OL potential, that will make it possible to reduce the model to the discrete 
form (see below), the latter condition amounts to the assumption that a z is essentially smaller that the lattice period. 
For instance, the trapping frequency cu z = 5 kHz yields a z ~ 0.5 /im for chromium atoms, while the OL period can be 
made equal to 5 /im, hence the latter condition can be readily satisfied. Actually, this means that the DD coupling may 
be considered as interaction between dipolar droplets with collective magnetic moments, trapped in local potential 
wells of the OL [H, [24| , while inside the strongly confined droplets the DD interaction effectively reduces to a contact 
form [H 

Continuing the derivation of the 2D model, axial width 7/, introduced above in Eq. is determined, in the lowest 
approximation by a local relation which reduces to an algebraic equation, 



rf = l + ^\(j>fr] 
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In the low-density limit, \4>\ « 1, Eq. (J6]) is tantamount to the GPE with the ordinary cubic local nonlinearity: 
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In this limit, Eq. ((?]) amounts to 77 pa 1. 

The choice of potential W (x, y) in the form corresponding to a deep OL suggests to approximate the wave function 
by a superposition of localized Wannier modes, $ m , n (x, y) 0, 4>(x, y) = Y^ m n 4>m,n&m,n (x, y). Projecting, as usual 
Q , equation ([6]) onto the set of this modes, Eq. ([6]) is reduced to the set of equations 
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where (m, n) play the role of 2D discrete coordinates replacing (x,y), the summation in the DD part excludes the 
term with (m',n') = (m,n), and C is an effective lattice coupling constant, determined by the overlap integral of 
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functions <E> m . n (a;, y) centered at two adjacent sites of the lattice. In addition to Eq. ([9]), discrete width function r) m . n 
obeys the respective counterpart of algebraic equation ([7]) , 

Vm,n = l + ^l0m,n| T)tn,n- (10) 

On the other hand, the 2D DGPE, which represents the discretization of Eq. ©, amounts to a single equation Q, 
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Stationary solutions to Eqs. ^ and (jlip with chemical potential fi are sought for as <j> m>n = u m n exp {—ijit), with 
discrete function w m .„ satisfying, respectively, the following stationary equations: 
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In the former case, r\ m ,n are related to u mi „ by Eq. (jTUJ> , with |0 m , n | 2 replaced by |u mj „| 2 . 

Generally, the expression for angles m . n in the term which describes the long-range DD interaction in the above 
discrete equations is cumbersome. However, following Refs. [H, [H, [Hj], it makes sense to focus on two most 
fundamental cases, with magnetic dipoles oriented either parallel or perpendicular to the (x, y) plane. In the former 
case, we choose the orientation of the dipole moments along the x. For the parallel orientation, the DD term in Eqs. 
(j!2[) and (|13[) is anisotropic, being attractive in one in-plane direction and repulsive in the other [l8| : 
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For the perpendicular orientation of the dipoles, the DD term is isotropic |15lll( 
expression: 
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taking a simple form, cf. the above 
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The character of the DD interaction in both cases is defined by the sign of T, which is related to the sign of coefficient 
g in the underlying GPE {T]). Normally, g is positive. However (as mentioned above), its sign may be reversed by 
means of a rapidly rotating magnetic field [13], allowing T to take negative values [HI, EH]. Thus, T > implies the 
repulsive isotropic DD interaction, while the anisotropic DD interaction is repulsive along discrete coordinate n, and 
attractive, being twice as strong in the absolute value, along to. For T < 0, the isotropic DD interaction is attractive, 
while in the anisotropic configuration the DD interaction is attractive along n and repulsive in the direction of to. 
The integral quantity that characterizes localized discrete modes is their norm (or power, in terms of optical models), 
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which is the dynamical invariant of Eqs. (J9} and (fTT]) . 
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For both in-plane and out-of-plane (perpendicular) dipole orientations, stationary equations (fT2|) and (fl3|) were 
solved by means of a numerical algorithm based on a modified Powell minimization method, which uses a finite- 
difference expression for the underlying Jacobian. In the analysis reported below, we focus on the three most funda- 
mental unstaggered localized modes ("staggered" are configurations with opposite signs of the field at adjacent sites, 
which is relevant to local models with the repulsive nonlinearity) : on-site modes, i.e., those centered on a lattice site, 
inter-site ones, which are centered between lattice sites in both directions m and n, and hybrid modes, centered on-site 
along one lattice direction, and between sites along the other. This nomenclature was adopted in previous studies of 
local discrete models Q. Results reported below were obtained using a square lattice composed of 11x11 or 21x21 
elements for on-site modes, and of 10x10 or 20x20 elements for inter-site and hybrid ones. It was checked that the 
results did not alter if an essentially larger lattices were used. 

III. FUNDAMENTAL UNSTAGGERED SOLITONS IN THE CASE OF ATTRACTIVE CONTACT 

INTERACTIONS 

The nonlinear localization of matter waves occurs in gaps of the linear spectrum of the underlying system. In the 
case of the attractive local nonlinearity, fundamental solitons populate the semi- infinite gap @, [29[ , which corresponds 
to regions [i < and // < 1, in the GPE and NPSE models, respectively. Discrete models, derived in the tight-binding 
approximation, cannot correctly describe the entire linear spectrum of the underlying continuum [301 ] . However, they 
can capture localized modes existing inside the semi-infinite gap. We will consider fundamental localized modes in 
this case (their topologically modified counterparts, such as discrete vortices [3ljj , will be considered elsewhere). 

To outline the entire existence region of the fundamental solitons, we followed the known approach, identifying 
regions where continuous-wave (CW) solutions with a given chemical potential, fi, are subject to the modulational 
instability, MI (it is assumed that the MI splits unstable CW states into arrays of solitons) [23|-[25j]. Further details 
of the respective numerical procedure are reported elsewhere [32j. It is relevant to notice that the instability of 
uniform states in the dipolar gas is determined by two different factors: the local self-focusing, which is the usual 
driving mechanism for the MI at small wavenumbers of perturbations, and, in addition, the roton instability at finite 
wavenumbers, induced by the long-range DD interactions [33| . 

A. The effect of the DD interaction on fundamental solitons 

As mentioned above, in the local model with the attractive contact interactions (H = —1), CW solutions are 
modulationally unstable in regions n < 1 and [i < for the DNPSE and DGPE models, respectively. Note that these 
values correspond to the upper boundaries of the semi-infinite gaps in the corresponding continuous NPSE and GPE 
models. The presence of the nonlocal DD interaction extends the MI domain of the CW states, thus enabling the 
existence of localized modes beyond the above-mentioned limit values, /i = 1 and /i = 0. 

In the case of the attractive contact interactions, families of fundamental unstaggered bright solitons of the three 
above-mentioned types - on-site, inter-site, and hybrid (see Fig. [T]) - have been found, as expected, precisely in the 
domain of the parameter space featuring the MI of the CW states. To this end, stationary equations (p~2|) and (Tl3|) 
were solved for the DNPSE and DGPE models, respectively. 

In the case of the attractive contact interaction, the attractive isotropic DD nonlinearity (that with the inverted 
sign) acts in the same direction and cannot change qualitative properties of the solitons. Therefore, only the case of the 
attractive on-site nonlinearity competing with its rep ulsive DD isotropic counterpart is really interesting, similar to 
the situation ID counterpart of the present model [20j , and in the Salerno model (a combined Ablowitz-Ladik/discrete 
nonlinear-Schodinger system) with competition between the attractive on-site and repulsive inter-site nonlinearities 
[34j . For this reason, in the case of the attractive contact interaction, we fix T > (the natural sign of the DD forces). 
In parallel, the interplay of the attractive on-site and anisotropic DD interactions will also be considered for T > 0. 

Basic characteristics of soliton families are presented by the P(/J,) dependence, which shows the norm versus the 
chemical potential. First, in Fig. Owe display P (fi) curves generated by the DNPSE (a) and DGPE (b) models for 
the on-site, hybrid and inter-site discrete solitons in the presence of the attractive contact interaction, while the DD 
interaction is still absent. These plots may be used as the reference point for the comparison with results obtained in 
the presence of the DD interactions. 

The P(p) dependencies for the soliton families affected by the isotropic (repulsive) and anisotropic DD interactions 
are displayed in Figs. [3]and[31 respectively. The comparison of the P(n) dependencies displayed in Figs. [3]and[5]with 
those in Fig. [2] shows that the repulsive isotropic DD interaction slightly extends the region in which the fundamental 
localized structures can be found, which is consistent with the results of the MI analysis for the CW solutions [32l |. 
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Simultaneously, the obtained P (fi) dependencies show that there is no qualitative difference between the DGPE and 
DNPSE models in this case. 

Usually, the Vakhitov-Kolokolov (VK) criterion, alias the slope condition, dP/dfi < 0, is considered as a necessary 
condition for the stability of soliton (35|, H|| . Although it is strictly applicable only to systems with a local power- 
law nonlinearity, it may be valid in more general situations too. Here we consider the slope criterion in both cases 
combining the cubic or nonpolynomial local nonlinearity with the nonlocal DD-interaction term. The slope criterion 
is violated at values of fi close to the right edge of the existence region of the fundamental solitons, see Figs. [H [3J 
andQ] 

The general condition which must be fulfilled for the stability of solitons is the absence of eigenvalues with positive 
real parts in the spectrum of small perturbations around the solitons (the spectral criterion); in fact, this means the 
eigenvalues must have zero real parts [361 ] . This criterion was implemented numerically by computing the corresponding 
eigenvalues, using linearized equations for small perturbations, similar to how it was done in Ref. Q. In the case of 
the on-site fundamental solitons, the real part of the eigenvalues vanishes in certain regions of parameter plane (C, (i), 
as shown in Figs. [5] (a) and (c) for both discrete models. This indicates the presence of a wide stability region for the 
on-site modes, as reported earlier in papers dealing with 2D solitons in the discrete nonlinear Schrodinger equation, 
without DD interactions Q, [I?! and (38|. The VK criterion is also satisfied in these cases. 

In the presence of the isotropic repulsive DD interaction the region of the existence of fundamental solitons expands. 
These extended regions are inside the domains where CW solutions are subject to the MI, in accordance with the 
above statement. Two significant effects are generated by the nonlocal DD interaction in this case. First, the isotropic 
repulsive DD interaction makes the stability regions wider for on-site fundamental solitons with respect to the model 
without the nonlocal interaction, as seen in Figs[5](b) and (d) for the DGPE and DNPSE models, respectively. Second, 
the nonlocal interaction generates new localized structures rising above a finite background. Properties of these new 
stationary structures, which are on-site centered, are considered in the following subsection. Here we only mention 
that the linear stability analysis shows that these solitons-on-the-background are unstable. On the other hand, in 
the presence of the anisotropic DD interaction, the stability domain of the on-site solitons shrinks (not shown here in 
detail). 

In the absence of the nonlocal interactions, a known result is that the fundamental solitons of the inter-site and 
hybrid types, which are characterized by greater values of the norms in comparison with the on-site modes, are 
unstable in the whole existence region Q ■ We have checked that DD interaction of either kind (isotropic repulsive or 
anisotropic) do not stabilize these modes. 

B. Solitons on a finite background 

As briefly mentioned above, a new feature revealed by the analysis in the presence of the repulsive isotropic DD 
interactions (i.e., those with the natural sign) is the role of the background in the formation of localized modes. 
This effect is most significant near the edge of the soliton's existence region (at /x close to or 1 in the DGPE 
and DNPSE models, respectively), where, in the presence of the isotropic DD repulsion, narrow localized modes are 
found, as stationary solutions, on top of a finite uniform background, see Fig. [HI We stress that such structures cannot 
be found in the absence of the DD interactions. When these interactions are present, the background may play a 
significant role in the formation of the localized structures, because the background as a whole couples to the central 
peak via the long-range DD forces. Actually, the long-range coupling gives rise to these newly formed solitons on the 
finite background (SFBs) when its strength exceeds a certain threshold level, see Fig. [7J The threshold becomes lower 
with the increase of coupling constant C, as seen in Fig. [7] Based on the waveform observed above the background, 
one can identify SFBs of the on-site, inter-site, and hybrid types. 

Localized SFB structures are generated by the MI of a finite-amplitude CW states in the nonlinear lattice with 
the DD interaction. In the course of the development of the MI, the interplay of the local (contact) and nonlocal 
(DD) nonlinearities arrests the exponential growth of the instability and imposes global correlations, which enable 
the creation of the SFB as an outcome of the nonlinear development of the MI. However, the SFBs themselves turn 
out to be unstable, eventually evolving into robust breathers, which consist of a localized vibrating peak and finite 
oscillating background, as shown below. 

C. Dynamical properties of the fundamental localized modes 

In addition to the modes described above, strongly pinned on-site solitons are found in certain parametric areas in 
the model with the attractive contact and isotropic repulsive DD interactions. Being very narrow, these modes are 
not significantly affected by the DD interaction. In general, these solutions correspond to early flat parts of the P(n) 
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plots in Figs. ©-((4]) (at fi < — 1 and fj, < in the DGPE and DNPSE model, respectively), and they have their 
counterparts ion the local model. These localized structures, such as the one displayed in Fig. ([H]) for the DNPSE, 
remain "frozen" under the action of arbitrary perturbations. 

The evolution of all the unstable fundamental solitons considered in the previous subsections follows approximately 
the same scenario, in the presence of the DD interactions. Direct simulations of Eqs. ([9]) and (|lip demonstrate that 
the outcome of the evolution under the action of perturbations is the emergence of persistent strongly pinned on-site 
breathers, surrounded by a finite oscillating background, see an example for the DNPSE in Fig. [5] The evolution may 
feature different transient stages, depending on the initial structure. For example, unstable on-site solitons directly 
evolve towards single-peaked tightly localized breathers, while inter-site solitons temporarily form localized structures 
with two peaks, see Fig. 1101 The eventually formed breathers are very robust and survive perturbations of any kind. 

A general conclusion is that the repulsive nonlocal DD interaction in the 2D discrete models does not prevent the 
formation of the extremely narrow localized structures, which are strongly pinned to the underlying lattice. These 
very narrow, dynamically stable modes have been associated with the quasi-collapse in the 2D lattice (3£^, [loT ]. 
Nevertheless, the DD interactions produce an essential effect, as they give rise to an oscillating background which 
surrounds the narrow peak, being coupled to it by the long-range DD interactions, see Figs. l9l and IT0| in the usual 
local model, with T = 0, the narrow peak would not be coupled to any background. 

With the increase the strength of the DD interaction, the evolution of the corresponding solitons of the SFB type 
becomes sensitive to the form of perturbations. Under small asymmetric perturbations (which actually correspond 
to small shifts of the initial soliton's location), or random perturbations, the SFBs decay. On the other hand, under 
small symmetric perturbations they evolve towards tightly bound single-peaked localized structures, surrounded by 
an oscillating background, similar to those displayed in Figs. l51and[TT)l 

The simulations demonstrate that the newly formed narrow breathers, built on top of the oscillating background, 
are robust modes which survive any type of perturbations. The intrinsic structure of these of the breathers depends 
on the initial perturbations which initiated the formation of the breather from an unstable soliton. Namely, for 
asymmetric or random perturbations, the repulsion exerted by the DD interactions is more significant than in the 
case of symmetric perturbations, making the emerging breather broader, with a richer internal structure, see Fig. 1111 

Finally, we have found that the fundamental solitons obtained in the model combining attractive contact and 
anisotropic DD interactions with a finite strength (recall they correspond to the in-plane orientation of the dipolc 
moments) are unstable in their entire existence region, except the on-site solitons in the region of very small C (the 
anti-continuum limit, as concerns the local coupling). We could conclude that the anisotropy of the DD interaction 
makes the solitons more sensitive to small perturbations (i.e., more unstable), in comparison to the case of the isotropic 
DD interactions. 



IV. FUNDAMENTAL SOLITONS GENERATED BY THE DD INTERACTIONS IN THE CASE OF THE 

REPULSIVE LOCAL NONLINEARITY 

In the local model with the repulsive contact interactions (H = +1), unstaggered localized modes cannot exist. 
However, the presence of the DD interactions, either attractive isotropic or anisotropic, opens a possibility to create 
such modes, as recently shown in discrete ID models 0, [25|, and was earlier predicted for quasi-2D continuous 
models in Refs. [15j (the competition of the repulsive contact and isotropic attractive DD interactions) and [18| (the 
interplay of the repulsive contact and anisotropic DD interactions). 

We have found that the long-range DD interactions between lattice sites may indeed create discrete 2D solitons in the 
case when the on-site nonlinearity is repulsive. Numerical calculations show that, as in the situations considered above, 
these discrete 2D solitons appear exactly in a part of the parameter space where the CW states are modulationally 
unstable (this is shown in detail in Ref. [32|]). This region belongs to the semi- infinite gap of the linear spectrum, which 
corresponds to [i < and /j, < 1, in the GPE and NPSE models, respectively. Again, three families of fundamental 
unstaggered discrete solitons - of the on-site, inter-site, and hybrid types - have been found in this situation as 
numerical solutions of stationary equations (fT2"|) and (JT3J) for the DNPSE and DGPE models, respectively. 

It is relevant to emphasize that the solitons are also obtained in the limit of the vanishing local interaction, i.e., the 
discrete solitons may be supported solely by the DD interaction, which is relevant to the experiment [ill ]. In other 
words, the attractive nonlocal DD nonlinearity is by itself sufficient for the formation of trapped localized modes in 
the 2D lattices. 



8 



A. Solitons generated by the attractive isotropic DD interaction 

In the case of the competition between the repulsive contact and attractive isotropic DD interactions, all unstaggered 
fundamental solitons have a bell-like shape, which is a consequence of the nonlocality of the interaction which creates 
them. This feature is also similar to what was found in the ID counterpart of the present models [23| . 

The P (n) curves for the on-site, hybrid and inter-site unstaggered solitons in the present case are much closer to 
each other than it was in the case of the attractive contact interaction, see Fig. [T2] On the other hand, numerical 
computations demonstrate a violation of the spectral stability criterion for all types of the unstaggered solitons in an 
almost entire existence region, except for the case of a small coupling constant, C < 0.1, and certain values of \i and 
r. With the increase of the strength of the isotropic DD attraction, the stability window shrinks with respect to the 
values of /i, as shown in Fig. [13] for C = 0.1. 

Direct numerical simulations show that the unstable bell-shaped localized structures evolve into the corresponding 
breathers, but, on contrary to the situations considered above, without significant generation of an oscillating back- 
ground, as shown in Fig. [14] (almost all the initial soliton's norm remains localized, rather than being transferred 
to the background, unlike the case of the attractive local interactions). The width of the localized breather remains 
nearly the same as that of the original unstable soliton. 

B. Solitons generated by the anisotropic DD interaction 

The anisotropic DD interaction can also provide for the formation of unstaggered fundamental discrete solitons 
of the on-site, hybrid and inter-site types in the case of the repulsive contact interactions. Examples of the P(p) 
characteristics for such discrete solitons at fixed T = 5 and two different values of the coupling constant, C = 0.02 
(close to the anti-continuum limit) and C = 2 (close to the continuum limit), are shown in Fig. 1151 The respective 
P(n) curves for solitons of the on-site and hybrid types are very close to each other. 

The linear-stability analysis indicates similar stability properties of the on-site and hybrid solitons. They are 
unstable in almost the whole existence region. As above, an exception is found close to the anti-continuum limit 
(C < 0.1), where stable solitons exist. The respective stability diagram is displayed in Fig. 02] (b). In this figure, the 
stable solitons of the on-site and hybrid types are found in the white window. The stability window for the inter-site 
solitons existing at C < 0.1 is smaller than for their on-site and hybrid counterparts. 

We have also considered the limit of large values of coupling constant C, i.e., the near-continuum limit, for which 
a narrow stability region for solitons had been reported in Ref. [l8l ]. produced by the variational analysis and 
direct simulations of the 3D continuous GPE. Taking into account the geometry, stable quasi-2D solitary modes, 
characterized by long-lived, slowly decaying oscillations were observed in certain parameter regions. However, in the 
discrete 2D DGPE and DNPSE models considered here (without the third dimension, that was explicitly present in 
the analysis presented in Ref. [H[), stable anisotropic solitons were not found in this part of the existence region. 

To summarize, in the model with the repulsive local interactions the long-range attractive isotropic or anisotropic 
DD interactions give rise to discrete solitons, but (marginally) stable ones can be found only close to the anti-continuum 
limit. This conclusion may be considered as a consequence of the long-range character of the DD interaction, which 
becomes more significant as the system is getting more discrete [27], [28| • 

V. CONCLUSIONS 

In this paper we have introduced the 2D discrete model for the dipolar BEC trapped in the deep optical lattice. We 
have analyzed the structure and dynamics of fundamental discrete solitons in two versions of the model, which are 
based, respectively, on the ordinary discrete form of the GPE (Gross-Pitaevskii equation) with the cubic nonlincarity, 
or the 2D discrete NPSE (nonpolynomial Schrodinger equation), both including nonlocal DD (dipole-dipole) interac- 
tions. Fundamental localized modes studied in this work are counterparts of solitons existing in the semi-infinite gap, 
in terms of the continuum limit of the discrete models. Only the situations with competing local and nonlocal inter- 
actions (attractive/repulsive, or vice versa) appear to be interesting. We have shown close similarity in the behavior 
of the fundamental unstaggered solitons in the discrete models of both the GPE and NPSE types in these cases. 

The presence of the repulsive isotropic DD interaction competing with the local attraction extends the existence and 
stability regions for fundamental on-site solitons. Due to the action of the nonlocal DD interactions, the background 
plays an essential role in the formation of a new type of unstaggered localized structures: near the edge of the 
semi-infinite gap of the corresponding continuous system, the generation of SFBs (solitons on a finite background) 
of the three types (on-site, inter-site and hybrid) was revealed by the systematic numerical analysis. These solitons 
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are unstable, but they spontaneously transform themselves into robust breathers, surrounded by a finite oscillating 
background. Actually, the presence of the background is a feature specific to the dipolar model. 

In the case of repulsive contact interactions, unstaggered fundamental solitons may exist due to the DD attraction. 
In the case of the isotropic attractive DD interaction, stable on-site, hybrid and inter-site solitons, with close values of 
the norm, can be found in the strongly discrete limit. In general, the fundamental solitons generated by the isotropic 
DD interaction are broad bell-shaped modes. Being unstable, they evolve into localized breathers with almost no loss 
of the norm. 
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FIG. 2: The dependencies for on-site, hybrid and inter-site solitons in the two-dimensional DNPSE (a) and DGPE (b) 

models, without the DD interactions; C = 2. 
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VI. FIGURES 
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FIG. 3: The -P(/x) dependencies for families of 2D solitons of the on-site (a), hybrid (b), and inter-site (c) types in the DNPSE 
model, in the presence of the isotropic DD interaction with relative strength F. The coupling constant is C = 2. 
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FIG. 4: The same as in Fig. [3] but in the case of the anisotropic DD interaction, which corresponds to the in-plane orientation 
of dipolar moments. 
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FIG. 5: (a,c) The stability diagram for on-site solitons in the (/i, C) plane in the case without the DD interaction: the DGPE 
(a) and NPSE (c) two-dimensional models. The influence of the isotropic repulsive DD interaction is presented in parts (b) and 
(d) for the DGPE and DNPSE models, respectively. The lattice coupling constant is C = 1. Solitons are stable and unstable 
in white and gray areas, respectively. 




FIG. 6: Comparison of an ordinary (fully localized) fundamental soliton (a) and SFB, i.e., the one sitting on a finite background 
(b) in the 2D DNPSE model. The powers (norms) of both solitons are the same, P = 3.5. The strength of the DD interaction 
is r = in case (a) and r = 1 in case (b), while the lattice coupling constant is the same, C = 2. The norm of the SFB, which 
includes the nonvanishing background, is finite due to the finiteness of the lattice. 
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Solitons on 
finite background 




FIG. 7: The existence diagram of SFBs (solitons on the finite background) in parameter plane (C, V) for the combination of 
the attractive contact and repulsive DD interactions in the 2D DNPSE model. The SFB is registered if the amplitude of the 
background exceeds 1/100 of the amplitude at the center of the soliton. With the increase of coupling constant C, the SFBs 
appear at a smaller strength of the DD interaction, F. 




FIG. 8: An example of a strongly pinned on-site soliton, which remains "frozen" (completely insensitive to perturbations) in 
the course of the evolution. Parameters are C — 2, fi — — 7, F = 1, which implies the combination of the attractive contact and 
isotropic repulsive DD interactions. In this and all following figures the results are presented for the nonpolynomial model. 
The findings are qualitatively the same in the model with the cubic nonlinearity. 




FIG. 9: The result of the application of a small perturbation to an unstable on-site soliton, which exists in the presence of the 
attractive contact and repulsive isotropic DD interactions (C = 2, /i = 0.89, F = 0.1). The soliton evolves into a persistent 
tightly bound structure, surrounded by an oscillating background. The snapshots correspond to t = 1 (a), t = 25 (b), and 
t = 50(c). 
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FIG. 10: The evolution of a perturbed inter-site soliton, for C = 0.2, \i — 0.7, V — 0.1, i.e., the combination of the contact 
attraction and isotropic DD repulsion. The emerging structure, as seen in panel (c), is a persistent one. The snapshots pertain 
to t = 1 (a), t = 25 (b), and t = 50 (c). 




FIG. 11: The evolution of an unstable perturbed soliton sitting on a finite background, for C = 2,[i = 1.1, F = 1, i.e., the 
combination of the attractive contact and isotropic repulsive DD interactions. 




FIG. 12: The P(n) dependencies, for C = 0.8 and T = —10, in the model with the repulsive contact interaction and isotropic 
DD attraction. For all the three types of the solitons - on-site, hybrid and inter-site - the P(/i) curves are very close to each 
other. 



16 




FIG. 14: The evolution of a perturbed bell-shaped on-site soliton with C = 0.8, fi = 0.5, T = —10, which corresponds to the 
combination of the on-site repulsion and isotropic DD attraction. Plots present snapshots at t — 1 (a), t = 10 (b), and t — 60 
(c). 




FIG. 15: The P(/J,) dependencies, for C = 0.02 (a) and C = 2 (b), with fixed T — 5, in the model with the repulsive contact 
and anisotropic DD interactions. In both plots, the P((i) curves for the three species of the discrete solitons are shown: on-site 
(solid line), hybrid (dashed line) and inter-site (dotted line). 



